Scalar radiation at finite temperature near a reflecting wall 
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This paper addresses the interplay between vacuum and thermal averages for massless scalar 
radiation near a plane wall of a large cavity where the Dirichlet boundary condition is assumed 
to hold. It is shown that the "Stefan constant" depends on the curvature coupling parameter and 
that the latter is restricted to a certain range by thermodynamic arguments. In more than three 
spacetime dimensions such a range contains the conformal coupling, but it does not contain the 
minimal coupling. Findings may be relevant in settings where arbitrarily coupled scalar fields are 
present. 

PACS numbers: 04.62,+v, 05.70.-a 



INTRODUCTION 

When physics of black body radiation in a cavity is de- 
rived in statistical mechanics textbooks the energy of the 
vacuum is usually ignored. The argument is that since 
the vacuum is by definition the state of minimum energy 
of a system, and one is interested only in the energy that 
can be extracted from the cavity, then the energy of the 
vacuum is effectively zero. In this context, electromag- 
netic radiation in thermodynamic equilibrium with the 
walls of a large cavity at temperature T is addressed ei- 
ther by considering it as an ideal gas of massless bosons 
or an ensemble of harmonic oscillators. Any of these tools 
leads to the energy density u and pressure P (equation 
of state) given by, 
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where the fundamental constants h and c in Eq. (Q]) point 
out the quantum-relativistic nature of the phenomenon. 
In another different but presumably equivalent ap- 
proach, one considers an ensemble of infinite cubic cavi- 
ties containing each the quantum field — quantum field 
at finite temperature. When techniques are applied to 
the electromagnetic field it results the following ensem- 
ble average for the stress-energy-momentum tensor, 
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with u and P given precisely by Eq. ([T}. If one of the 
walls of the cavity is brought isothermally from infinity 
and parallel to the opposite fixed wall, (T^ v ) changes rad- 
ically, loosing its isotropic character although remaining 
uniform [l|. When the distance between these two walls 
is small enough, (T^ w ) is such that hardly depends on T; 
in fact it does not even vanish when T — 0. At this stage, 
the two walls are attracted to each other — the Casimir 
effect — and that is a surprising dynamical manifestation 



of the vacuum state in quantum field theory (see reviews 
01). 

There are also manifestations of the electromagnetic 
vacuum when the walls are far apart. When a large cav- 
ity has curved walls [J-tZ| , the behavior of electromagnetic 
radiation in the bulk is still essentially given by Eqs. @ 
and (JXJ) . Approaching a wall one sees that ( i T fil> ), in fact, 
is not uniform. In particular near the wall the expression 
for u in Eq. ([T]) becomes a subleading contribution in 
(T 00 J. The corresponding leading contribution is a tem- 
perature independent term, which diverges as the curved 
wall is approached if theperfect conductor boundary con- 
ditions are taken on it [4J, [7| . 

Considering now massless scalar radiation, deep in the 
bulk of a large cavity Eqs. (0) and ^ hold, after halving 
u due to absence of polarization. In analogy with elec- 
tromagnetic radiation in a cavity with curved walls, one 
will eventually notice that (T^") is not uniform by mov- 
ing toward a wall, but this time regardless of whether 
the wall is curved or not. As has been shown in Ref. 
[6J, the thermal behavior of (T 00 ) near a reflecting plane 
wall will be still Planckian, but the corresponding "Stefan 
constant" will depend on the kind of boundary condition 
— Dirichlet or Neumann — that is assumed to hold on 
the wall. 

Another new ingredient that thermal scalar radiation 
brings about, and that has been overlooked in the liter- 
ature, is the dependence of (T Ml/ ) on the curvature cou- 
pling parameter £. As is well known, in flat spacetime 
£ does not appear in the wave equation; but it does ap- 
pear in the corresponding stress-energy-momentum ten- 
sor. The reason is that the variation of the action with re- 
spect to the metric that leads to T^ v is made before solv- 
ing the Einstein equations which yield the flat geometry. 
Using the scalar Feynman propagator to obtain (T^ v ) ac- 
cording to the "point splitting procedure" [8j , this paper 
investigates the dependence of (T^ v ) on £ near a reflect- 
ing plane wall of an infinite cavity in flat spacetime with 



arbitrary number of dimensions. The main result is that 
by requiring thermodynamic equilibrium only values of £ 
restricted to a certain range are acceptable. In the rest 
of the text, unless stated otherwise, Kb = ft = c = 1. 



THE SCALAR PROPAGATOR 

A plane wall is taken at x = 0, where the massless 
scalar field <j> vanishes — Dirichlet's boundary condition. 
All the other walls of the (N — l)-dimensional cavity are 
at infinity. Radiation in the cavity is assumed to be in 
thermodynamic equilibrium with the wall at temperature 
T, so one uses the formalism of analytically continuing 
time t to imaginary values, taking it periodic with period 
1/T. Observing these prescriptions the equation [8[ 
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is solved to obtain the scalar propagator. The solution 
is well known for N = 4 (see Ref. [8], and references 
therein). When N > 3, 
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where a± := (t - f - in/T) 2 - (x ± x') 2 - (y - y'f - 
(z — z 1 ) 2 — ■ ■ • plus an infinitesimal negative imaginary 
term. Noting that the term corresponding to n — and 
involving <r_ in Eq. Q is simply the zero temperature 
propagator in ordinary Minkowski spacetime, direct ap- 
plication of D x to Eq. yields promptly Eq. ([3]). In 
the following, the Minkowski propagator will be dropped 
in order to implement renormalization. It is worth men- 
tioning that for N < 3 familiar divergences arise [9(. As 
will be seen shortly, such divergences do not bother (T^") 
though when N > 2. 



THE STRESS-ENERGY-MOMENTUM TENSOR 

The expectation value of the stress-energy-momentum 
tensor can be obtained by applying the differential oper- 
ator m 

v p,v ._ (j _ 2£)d fi d v ' - 2S,d"d v + (2£ - l/2)ri^d x dx' 
to the renormalized propagator, 
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resulting that (T* 1 ") is conserved, traceless when £ = £jv 
(conformal coupling) , 
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and with diagonal form, 

(TH=diag(p ) Pj.,P||,--- ) P||). 
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The only component that is uniform in Eq. ([6]) is the 
pressure perpendicular to the wall, 
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The interest here is on the asymptotic behaviors of the 
quantities in Eq. ©. For convenience formulas will be 
specialized to non negative x. In the bulk (Tx 3> 1) one 
finds for the energy density, 



p=(N- 1)Pjl + Pc 
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and for the pressure parallel to the wall (N > 3), 
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where P c i aS s is given by the negative of Eq. ((9]) with 
£ - 1/4 replaced by £ - &v-i- Eqs. © and JTOJ) are 
exact up to exponentially small corrections, and at high 
temperatures they apply also near the wall. 

Moving now near the wall (Tx C 1), P± is still given 
by Eq. (|7|) since it is uniform. Neglecting terms with 
higher powers of Tx, it results that 
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Clearly, Eqs. (jll|l and (|12j) hold also everywhere at low 
temperatures. 

Part of the material in Ref [6] concerns with thermal 
scalar radiation near a reflecting plane wall. By setting 
in Eq. flS) N = 4 and £ = (minimal coupling), or £ = 
1/6 [conformal coupling when N = 4, cf. Eq. |[5J], the 
results in Ref. |6| are successfully reproduced. Eq. (fTTj) 
also reproduces the corresponding expression in Ref. [6j , 
when TV = 4 and £ = 1/6. Ref. [lOj has investigated 
thermal scalar radiation between two parallel reflecting 
walls separated by unity, providing a formula for (T M,y ) 
with arbitrary £ and N = 4. It has been checked that, by 
reintroducing arbitrary distance between the walls and 
taking it to infinity, the resulting asymptotic behaviors 
of (T^ v ) in Ref. [l0( agree with those above, after setting 
N = 4. 



DISCUSSION 



[cf. Eq. CP], 



The expression for P±_ in Eq. ([7]), which holds every- 
where at arbitrary T, is precisely the usual black body 
radiation pressure in TV dimensions. Deep in the bulk 
(x — > oo) Pciass and P c i aS s in Eqs. §E§ and (TTU)) can 
be neglected, resulting the usual relations proper of the 
uniform and isotropic black body radiation. 

By reintroducing dimensionful h in Eqs. ((5J) and (|10p. 
one sees that the "classical" corrections /5 c i aS s and P c iass 
carry H° whereas the black body quantities carry a neg- 
ative power of h, namely, h l ~ N . Although the termi- 
nology "classical" is at some extent justified, by setting 
h — >• the "classical" corrections effectively vanish since 
the black body quantities diverge. In the context of 
TV = 4, Ref. [6j has pointed out that the only relevant 
correction to the black body contribution in Eq. © is 
linear in T. Now it can be seen from Eq. ([9]) that this 
is so regardless the number of dimensions N. This fact 
resembles the equipartition principle of energy, and it 
could not be reached simply by using dimensional argu- 
ments. At this point it should be noted that "classical" 
contributions also appear in the context of Casimir's ef- 
fect at finite temperature [l|,|j| (see also Refs. [2j,|3|, and 
references therein). 

Whereas only the subleading contribution in the ex- 
pression for p in the bulk [cf. Eq. (jHJ)] carries dependence 
on £, both terms in Eq. (fTTj) depend on it. As the wall 
is approached, the vacuum energy density in Eq. (fTTj) 



diverges for £ ^ £jv 11]: and this limitation is commonly 



understood as a consequence of replacing a real reflect- 
ing plane wall (which has some thickness) by a boundary 
condition on a plane [5J (see also Ref. [6j , for an alterna- 
tive interpretation). Although the subleading contribu- 
tion in Eq. (fTTj) does not depend on x, strictly speaking it 
should not be taken as purely thermal. Its mixed nature 
(vacuum-thermal) is reflected in the fact that it depends 
both on £ and T. As has been pointed out in Ref. [6j for 
N — 4, and mentioned earlier in the text, the behavior 
in Eq. |[TT)) resembles that for electromagnetic radiation 
near a curved wall where the perfect conductor boundary 
condition is assumed to hold [4j . An important difference 
though is that the subleading contribution in the case of 
the electromagnetic radiation is the very black body en- 
ergy density, and thus purely thermal. Another feature 
near the wall (or everywhere at low temperatures) that 
should also be remarked is that whereas P± in Eq. (|7|) 
is always positive, each term of Pii in Eq. (|12l) can be 
negative depending on £. 

The expression for the specific heat per unit of volume 
in the bulk [cf. Eq. (|SJ)] is given by 
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where the factor multiplying NT N 1 is the familiar "Ste- 
fan constant" in N spacetime dimensions. Near the wall 



C 'V = ^ r (f ) CW(1 - ^)T N -\ (14) 

which also has Planckian form, but now the "Stefan con- 
stant" depends on £. Examining these equations one sees 
that cy in the bulk is always positive, but near the wall 
it can be negative. Assuming scalar radiation in ther- 
modynamic equilibrium on each side of the wall, if the 
temperature Tr on the right side is slightly higher than 
that on the left, Tr > Tl, the net pressure on the wall 
will be given by the difference of the corresponding val- 
ues of P± as in Eq. (|7J). If a tiny puncture is made 
on the wall, then energy will flow from the right to the 
left — these are consequences of the integral version of 
the conservation law T^ v ,„ — 0, and that the energy flux 
density is equal to the momentum density. Noting Eq. 
(fTT)) . when £ > 1/4, cy < resulting that Tr becomes 
even higher and Tj, lower, i.e., the system is carried way 
from thermodynamic equilibrium. It follows then that £ 
must not be greater than 1/4. 

A similar kind of arguments near the wall show that £ 
may also have a lower bound. One begins by considering 
a tiny cubic region of the scalar radiation near the wall, 
and say facing the wall. It is further assumed that the 
temperature Tj„ inside the cube is slightly higher than 
the temperature T out outside. Using once again the inte- 
gral version of the conservation law T^ v ,„ = and that 
the energy flux density is equal to the momentum den- 
sity, it follows from Eqs. ©, Q, dUJ and JJJ} that the 
flux of energy out of the tiny cubic region is given by 
[1 - (TV - 2)(1 - 4£)]{TF - T* t ) times a positive factor 
that includes a small time interval. As Ti n > T out and 
£ < 1/4 by assumption, the energy inside the cube must 
decrease [cf. Eq. ([14"])]. otherwise the system will run 
way from thermodynamic equilibrium. Accordingly, the 
flux above must be positive, i. e., 1 — (N — 2)(1 — 4£) > 0. 
Putting all together, £ needs to be such that 



for N = 2, and 
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for N > 2. It is rather curious that the lower bound in 
Eq. (1161) is the conformal coupling in one less dimension 
[cf. Eq. J5J], and that the upper bound is the limit of £at 
as N — > oo. Thus the conformal coupling £/v fits in Eq. 
(fTBl) . but the minimal coupling £ = does not. As can 
be readily seen, to eventually including equalities in Eqs. 
(fTS")) and (fll)|) , further corrections in Eqs. ([TT]) and ([12"]) 
are required. It should be additionally noted that when 
N — > oo, the interval in Eq. (|16p narrows to a point, 
encapsulating the conformal coupling. 



A connection between Eq. (|T6|) and the "classical" 
corrections in Eqs. © and (jlOp is worth pointing out. 
The "classical" corrections are related by 
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and therefore Eq. (|16|) implies that P c iass and /O c iass 
have the same sign. Now, whereas Eq. (IT7|) is a feature 
of the behavior of the scalar radiation in the bulk, the 
permissible range in Eq. (TT51) arises by considering its 
behavior near the wall. 

Mathematically speaking, a wall is a boundary con- 
dition assumed to hold on a surface to prevent flux of 
energy across the surface. As Eqs. (fT5*)) and (fTtj)) were 
obtained assuming the Dirichlet boundary condition, it is 
pertinent to ask if another commonly used boundary con- 
dition, namely the Neumann boundary condition, leads 
to the same result. In other words one may wish to know 
how much dependent on the type of boundary condition 
Eqs. (I15[) and (J16I) are. A preliminary answer can be 
obtained from Ref. [lOj . After manipulating formulas 
in Ref. [10| concerning the Neumann boundary condi- 
tion, one finds as permissible range —5/4 < £ < 7/8, 
which is to be compared with Eq. Q16p for N = 4, i.e., 
1/8 < £ < 1/4. It follows that, in this case, the per- 
missible range corresponding to the Neumann boundary 
condition not only includes the minimal coupling £ = 0, 
but also contains the range corresponding to the Dirichlet 
boundary condition. The Dirichlet boundary condition 
is therefore more restrictive. 

An important issue to be addressed here is to which 
extend the atmosphere near the wall can be considered 
as an ideal gas of massless scalar bosons at temperature 
T. Taking for convenience N = 4, one may ask to which 
extend a standard formula such as (reintroducing dimen- 
sionful constants) 
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still holds in the present situation [cf. Eq. (jHJ]. It should 
be noted that the right hand side of Eq. (fi"8|) is in fact 
equal to U/V, and strictly speaking the equality might 
hold only if the internal energy U were uniformly dis- 
tributed in the cavity, which is clearly not the case for the 
scalar radiation. Nevertheless, deep in the bulk (x — > oo) 
Eq. ([18]) and its corresponding equation of state p = p/3 
hold as follows from the discussion opening this section. 
Similarly one might guess that the temperature depen- 
dent term in Eq. (|lll) [see also Eq. (|14l) ]. which holds 
near the wall, could be associated to an effective Planck's 
distribution as in Eq. (|18|) . However, recalling how pres- 
sure of a gas on the walls of a cavity is obtained from 
the distribution function, one sees clearly that such an 
effective Planck's distribution would fail in reproducing 
the pressure in Eqs. (0) and (fi"2j) . This conclusion has 



to do with the already mentioned mixed nature (vacuum- 
thermal) of the scalar radiation at finite temperature near 
the reflecting wall. 

Perhaps this duality can be better appreciated by say- 
ing, as in Ref. [4], that the atmosphere near the wall is 
a mixture of virtual and real bosons. Indeed, if initially 
at any point in the cavity the specific heat is given by 
the nearly black body expression in Eq. (I13[) . and then 
the temperature is lowered such that eventually cy will 
be given by Eq. (fill), as T -► 0, cy ->• and P± -)■ 0, 
resulting that at T = only virtual bosons are left in the 
cavity from where no further energy can be extracted — 
only vacuum energy is left behind. 

Differently from the scalar radiation, the thermal be- 
havior of the electromagnetic energy density p in the bulk 
and near a curved wall of a large cavity is the same [4| . 
Thus an expression like (precisely twice that in) Eq. (|T5|> 
yields the thermal contribution in the bulk and also near 
the wall, corresponding to the usual radiation pressure 
[cf. Eq. (JXJ)] . It follows then that, in this case, one can 
think of an ideal gas of massless bosons (photon gas) in 
the bulk, as well as near the wall. 

Before closing, it should be noticed that by dropping 
the vacuum energy density, one can check that the inte- 
gral over space of p in Eq. © does not depend on £, for 
N > 2. This has been pointed out also in Ref. [lfj when 
N = 4, for a configuration of two parallel walls, and it 
is consistent with the fact that standard black body ex- 
pressions do not carry dependence on £. 

Recapitulating, this paper investigated the vacuum- 
thermal nature of the scalar radiation. The main con- 
clusion is that statistical mechanics of a scalar field re- 
stricts the values of the curvature coupling parameter £. 
Corresponding to the Dirichlet boundary condition on a 
plane wall, such a restriction is characterized by the in- 
equalities in Eqs. (fl5)) and (p~6|) . Further still unclear 
role seems to be played by the bounds in Eq. (p~6|) . For 
example, the symmetrical appearance of these bounds in 
the "classical" corrections to the black body expressions 
in the bulk asks for more investigation. 
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